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Abstract. In the recent years we have performed various experiments on the collision dynamics of highly porous dust 
aggregates and although we now have a comprehensive picture of the micromechanics of those aggregates, the macroscopic 
understanding is still lacking. We are therefore developing a mechanical model to describe dust aggregate collisions with 
macroscopic parameters like tensile strength, compressive strength and shear strength. For one well defined dust sample 
material, the tensile and compressive strength were measured in a static experiment and implemented in a Smoothed Particle 
Hydrodynamics (SPH) code. A laboratory experiment was designed to compare the laboratory results with the results of the 
SPH simulation. In this experiment, a mm-sized glass bead is dropped into a cm-sized dust aggregate with the previously 
measured strength parameters. We determine the deceleration of the glass bead by high-speed imaging and the compression 
of the dust aggregate by x-ray micro-tomography. The measured penetration depth, stopping time and compaction under the 
glass bead are utilized to calibrate and test the SPH code. We find that the statically measured compressive strength curve is 
only applicable if we adjust it to the dynamic situation with a "softness" parameter. After determining this parameter, the SPH 
code is capable of reproducing experimental results, which have not been used for the calibration before. 
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INTRODUCTION 

Planets form from micrometer-sized dust grains, collid- 
ing at low velocities and sticking to each other due to at- 
tracting van-der-Waals forces [1]. From the interplay be- 
tween laboratory experiments [2] and molecular dynamic 
simulations [3, 4] we have a good picture of the micro- 
physics of those processes, but as aggregates grow larger 
(e.g. 1 mm), the understanding of the collision mechanics 
is severely lacking. For these sizes (and the correspond- 
ing velocities), collisions between equal sized aggregates 
rather lead to bouncing or fragmentation than to sticking. 

Millimeter-sized aggregates - as precursors of plan- 
ets - are expected to be highly porous, having a volume 
filling factor (the volume fraction of material) of only 
few percent up to few ten percent. This strongly cohesive 
material is comparable to millimeter-sized dust clumps 
present in conventional dry powders (e.g. cocoa, confec- 
tioners' sugar) but probably being more porous. One well 
analyzed analog material are dust aggregates formed by 
the random ballistic deposition (RBD) method, experi- 
mentally introduced in Ref. [5]. These dust aggregates 
have a volume filling factor of = 0. 15 (for monospheric 
dust) and are produced in our laboratory in macroscopic 
2.5 cm samples. Using these dust samples, we performed 
various collision experiments to study the further growth 
and evolution of planetary bodies [e.g. 1, 6]. However, 
as experiments cannot be performed for any relevant set 



of parameters (e.g. collisions of meter-sized bodies), a 
collision model is required to cover the wide parameter 
space occuiTing for protoplanetary dust-aggregate colli- 
sions, i.e. dust aggregate sizes of up to 1 km and col- 
lision velocities in the range of 10"^ - 10^ ms-'. The 
approach is therefore to measure macroscopic parame- 
ters of the laboratory dust samples and implement them 
in a numeric simulation of a dust aggregate collision, 
which is concurrently performed in a laboratory exper- 
iment. Comparing experiment and simulation and defin- 
ing and determining the free parameters, yields a cali- 
brated collision model to perform collision simulations 
with parameters unaccessible to laboratory experiments. 

SMOOTHED PARTICLE 
HYDRODYNAMICS FOR DUST 
COLLISIONS 

For the simulation of dust aggregate collisions we use the 
Smoothed Particle Hydrodynamics (SPH) method with 
extensions for the treatment of solid and porous media. 
A comprehensive description of the meshless Lagrangian 
particle method SPH can for example be found in [7]. In 
this scheme the continuous solid objects are discretized 
into interacting mass packages ("particles") caiTying all 
relevant continuous quantities. Time evolution is com- 



puted according to the Lagrangian equations of contin- 



uum mechanics: 
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Here, p denotes the density, v the velocity, D the dimen- 
sion and Gafi the stress tensor, defined as 
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It consists of a pressure part with pressure p and and a 
shear part given by the traceless deviatoric stress tensor 
Safj. Its time evolution is modelled according to Ref. [8]. 
This set of equations is closed by a suitable equation of 
state and describes the elastic behavior of a solid body. 
Together with a suitable damage model, which we do not 
adopt, the authors of Ref. [8] have modelled collisions 
between brittle basaltic rocks using this scheme. 

In contrast to this, we simulate the plastic behavior of 
porous bodies. Therefore we adopt a modified version 
of the porosity model by Sirono [9]. According to this 
approach, plasticity is modelled within the equation of 
state and porosity is given by 1 — p /po, where p denotes 
the actual and po the bulk density of the material. The 
pressure is limited by the compressive strength E(p) as 
upper bound and the (negative) tensile strength T{p) 
as lower bound. In between, the solid body experiences 
elastic deformation, whereas outside this regime it is 
deformed plastically. Thus, the full equation of state 
reads 



p{p)={ Kip'a)ip/p'o-l) 
T{p) 



P>Pc 
- < P < Pe+ 
P <Pc 
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The quantity p'q denotes the density of the material 
at zero external stress. p+ and p^T are limiting quanti- 
ties, where the transition between the elastic and plastic 
regime for compression and tension, respectively, takes 
place. Once a limit is exceeded, the material leaves the 
elastic path where energy is conserved, and loses inter- 
nal energy by following the paths of the compressive and 
tensile strength. 

In a previous work, Sirono [9] adopted power laws 
from measurements with toner particles for E(p), T{p) 
and the bulk modulus K{p ) in order to simulate porous 
ice, which is a crude extrapolation. For our approach, we 
used the material properties, measured for well defined 
RBD dust aggregates [5]. For the tensile strength we 
adopted the measurement of Ref. [5] who measured the 
tensile strength for porous ((j) =0.15) and compact = 
0.54) aggregates and found an agreement with a linear 
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FIGURE 1. Compressive strength curve for unidirectional 
(dashed line), omnidirectional (solid line), and dynamic (dotted 
line) compression. The inset illustrates the setup for the omni- 
directional compression measurement. 



dependence between tensile strength and numbers of 
contact per cross-sectional area. This yields a tensile 
strength of 



r((/)) = -(io 



2.8+1.. 



Pa. 



(5) 



For the compressive strength we started with the com- 
pression curve 0(2) measured in Ref. [5] (unidirectional, 
ID compression) and also made a new compression mea- 
surement (omnidirectional, 3D, see inset in Fig. 1, [10]). 
Both compressive strength curves are displayed in Fig. 1 . 
They resemble in shape and can be described by 
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with four free parameters 0i , 02, A, and pm- For unidirec- 
tional compression we found pm = 5.6 kPa, 0i = 0.15, 
02 = 0.33, and A = 0.33, while for ominidirectional 
compression measurements performed in this work the 
parameters are identified as p^ ~ 13 kPa, 0i — 0.12, 
02 = 0.58, and A = 0.58. For unidirectional, static com- 
pression, the material can creep sideways, releasing pres- 
sure. As we did not observe this in the dynamic experi- 
ments (next section. Fig. 2), we expect omnidirectional 
compression, defining the parameters 0i and 02. For low 
pressures, A • In 10 is the slope of a power law found by 
divers authors [5, 11] and does not leave very much mar- 
gin. Thus, we take as the only free parameter with 
the most influence on material softness, shifting the com- 
pression curve towards lower pressures and determine 
Pyy, within the calibration procedure. As the the shear 
strength was not measured so far, we follow Sirono [9] 
and take Y = y/L\T\. 
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FIGURE 2. The density plot reveals the compaction of the 
dust sample under the glass bead (red, saturated) with a volume 
filling factor of 0.20 to 0.25 (green) in a well confined volume 
of approximately one sphere volume. The original dust material 
has a volume filling factor of 0. 15 (light blue). 

CALIBRATION EXPERIMENT 

In the experiments, a solid projectile was dropped into 
a 2.5 cm diameter, highly porous dust aggregate con- 
sisting of 1.5 fim Si02 monomers (RBD dust aggregate 
from [5]). The experiments were performed in vacuum 
(0.1 mbar) such that gas effects did not play a role. The 
projectile was either a glass bead of 1 mm diameter or a 
cylindrical plastic tube with a 1 or 3 mm diameter epoxy 
droplet at the bottom (representing a 1 or 3 mm glass 
bead). The epoxy projectile had a mass corresponding 
to a glass bead and therefore was longer and could be 
observed for a penetration deeper than its diameter. For 
15 experiments with elongated projectiles, the decelera- 
tion curve and, thus, the penetration depth and the stop- 
ping time were measured with a high-speed camera. The 
stopping time for one projectile diameter was found to be 
rather constant, (3.0 ±0.1) ms for 1 mm projectiles and 
(6.2 ±0.1) ms for 3 mm projectiles, while the penetra- 
tion depth depended on the projectile size and the impact 
velocity (drop height). Details on the full deceleration 
curve can be found in [10]. The penetration depth could 
well be approximated by 

/ . m^ s \ mv 

where v, m, and A — nR^ are the projectile velocity 
at the time of first contact, the projectile mass, and its 
maximum cross-sectional area (see also Fig. 4 in the 
next section). The stopping time and the relation for the 
penetration depth will be used for the calibration later on. 

Two experiments with a 1 mm glass bead impact- 
ing with (0.8 ±0.1) m s^' were analyzed using an x- 
ray micro-tomograph (Micro-CT SkyScan 1074). In this 
method, the dust sample with the embedded glass bead 
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FIGURE 3. Varying the parameter in the simulation 
changes the softness of the material and yields different pen- 
etration depths and stopping times. The best agreement with 
the experiments was found for pm = 1.3 kPa. 



was positioned between an x-ray source and a detector 
and rotated stepwise. Based on the 400 resulting trans- 
mission images, a 3-dimensional density reconstruction 
was calculated using the SkyScan Cone-Beam Recon- 
struction Software. Assuming cylindrical symmetry in 
the axis of penetration, the density was averaged to one 
vertical section and displayed in the contour plot in Fig- 
ure 2. The compaction of the dust is clearly visible in a 
confined volume under the sphere. The green color mark- 
ing this volume denotes a volume filling factor of 0.20 to 
0.25, while the material around this compaction zone is 
virtually unaffected with an original volume filling factor 
of ^ w 0.15. The distribution of compacted material will 
be used for the verification of the SPH code. 



CALIBRATION AND VERIFICATION OF 
THE SPH CODE 

In the 2D SPH simulation, an infinite cylinder with 
1 . 1 mm diameter (p = 2540 kg m^^*) impacts into a dust 
sample with a cross section of 8 x 5 mm^ at a velocity 
of 0.65 m s^'. For comparison with the 3D experiments 
a correction factor of for the penetration depth is re- 
quired (see [10]). For the calibration we will use the pen- 
etration depth, which has to be 0.82 mm (Eq. 7 and cor- 
rection factor), the stopping time of 3 ms, and the com- 
paction under the glass bead from the x-ray experiments. 

Using the shear strength formalism introduced by 
Sirono [9], e.g. Y = y/L\T\, we have one free parame- 
ter with influence on the outcome of the simulations, i.e. 
the adjustment of the compressive strength curve for dy- 
namic compression by the parameter (see Eq. 6). A 
detailed study on the shear model can be found in [10]. 




FIGURE 4. The penetration depth as a function of momen- 
tum per cross-sectional area can well be reproduced with the 
SPH code. 



The parameter pm defines the softness of the material 
and decreasing this parameter yields deeper penetrations. 
Figure 3 describes the calibration of this parameter. The 
horizontal line denotes the expected penetration depth 
from Eq. 7, while the vertical line represents a mean 
stopping time of 3 ms for 1 mm spheres. A value of = 
1.3 kPa yields the best agreement with the experiments 
and thus we will use this value for further tests. 

One test for the validation of the SPH code is the 
penetration depth relation for different sizes and different 
velocities, reproducing Eq. 7. The penetration depths 
of the experiments and the penetration depths of the 
simulations (with correction factor) are plotted in Fig. 4. 
For mvA^^ ?i 1 kg m^' s^' experiment and simulation 
are in very good agreement and the SPH code succeeds 
in the scaUng of radius and velocity. A second vaUdation 
test is the comparison of the compressed volume with the 
results of the x-ray measurement (Fig. 2). From the 3- 
dimensional density data, one can determine the volume 
fraction that is compressed to a given volume filling 
factor, which is plotted in Fig. 5. The values for (j) < 0.2 
represent the original material being distributed around 
the mean filUng factor of 0.15. Also for this test, we find 
a good agreement between experiment and simulation. 

CONCLUSION 

We have developed an SPH code for the description of 
colUsional interaction between high-porosity dust aggre- 
gates. The code was calibrated by static and dynamic lab- 
oratory experiments using macroscopic RBD aggregates 
with an uncompressed filling factor of ^ = 0.15. The cal- 
ibrated SPH code correctly predicts the size and veloc- 
ity dependence of the penetration depth for an impacting 
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FIGURE 5. Experimental results (symbols) and simulation 
(solid line) of the compressed volume. 



solid projectile as well as the compressed volume. 



ACKNOWLEDGMENTS 

We thank M.-B. Kallenrode and the University of Os- 
nabriick for providing access to the XRT setup. Simula- 
tions were performed on clusters of the computing cen- 
ter (ZDV) of the University of Tiibingen. This project is 
funded by the Deutsche Forschungsgemeinschaft within 
the Forschergruppe 759 "The Formation of Planets: The 
Critical First Growth Phase" under grants Bl 298/7-1, Bl 
298/8-1, and Kl 650/8-1. 



REFERENCES 

1. J. Blum, and G. Wurm, Annual Review of Astronomy and 
Astrophysics 46, 21-56 (2008). 

2. J. Blum, and G. Wurm, Icarus 143, 138-146 (2000). 

3. D. Paszun, and C. Domiruk, Astronomy and Astrophysics 
484, 859-868 (2008), 0802. 1832. 

4. C. Dominik, and A. G. G. M. Tielens, The Astrophysical 
Journal 480, 64^ (1991). 

5. J. Blum, and R. Schrapler, Physical Review Letters 93, 
115503 (2004). 

6. D. Langkowski, J. Teiser, and J. Blum, TTie Astrophysical 
Journal 675, 764-776 (2008). 

7. J. J. Monaghan, Reports of Progress in Physics 68, 
1703-1759 (2005). 

8. W. Benz, and E. Asphaug, Icarus 107, 98 (1994). 

9. S.-I. Sirono, Icarus 167, 431^52 (2004). 

10. C. Guttler, M. Krause, R. Geretshauser, R. Speith, 
and J. Blum, The Physics of Protoplanetesimal Dust 
Agglomerates. IV. Towards a Dynamical Collision Model 
(2009), submitted. 

11. J. M. Valverde, M. A. S. Quintanilla, and A. Castellanos, 
Physical Review Letters 92, 258303 (2004). 



